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ABSTRACT 

We show how wavelength-dependent time delays between continuum flux varia- 
tions of AGN can be used to test the standard black hole-accretion disc paradigm, 
by measuring the temperature structure T(R) of the gaseous material surrounding 
the purported black hole. Reprocessing of high energy radiation in a steady-state 
blackbody accretion disc with T oc i? -3 / 4 incurs a wavelength-dependent light travel 
time delay r oc A 4 / 3 . The International AGN Watch multiwavelength monitoring cam- 
paign on NGC 7469 showed optical continuum variations lagging behind those in the 
UV by about 1 day at 4800A and about 2 days at 7500A. These UV/optical con- 
tinuum lags imply a radial temperature structure T oc i? -3 / 4 , consistent with the 
classical accretion disc model, and hence strongly supports the existence of a disc 
in this system. We assume that the observed time delays are indeed due to a clas- 
sical accretion-disc structure, and derive a redshift independent luminosity distance 
to NGC 7469. The luminosity distance allows us to estimate a Hubble constant of 
Hq (cosi/0.7) 1 ^ 2 = 42 ± 9kms -1 Mpc -1 . The interpretation of the observed time 
delays and spectral energy distribution in the context of an accretion disc structure 
requires further validation. At the same time, efforts to minimize the systematic uncer- 
tainties in our method to derive a more accurate measurement of Hq, e.g by obtaining 
an independent accurate determination of the disc inclination i or statistical average of 
a moderate sample of active galaxies, are required. However, this remains a promising 
new method of determining redshift-independent distances to AGNs. 
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1 INTRODUCTION 

The Hubble constant Ho and deceleration parameter go are 
fundamental parameters in standard cosmology, measuring 
respectively the rate at which the Universe is expanding and 
the rate at which that expansion is impeded by the attrac- 
tive force of gravity. Moreover, Ho determines a size scale 
and age of the Universe, constrains the baryonic density pro- 
duced in the Big Bang, the amount of dark matter in the 
Universe, and the epoch for galaxy and quasar formation 
in the early Universe. A measurement of the deceleration 
parameter constrains the geometry of the Universe. 

The value of the Hubble constant remains in dispute af- 
ter over half a century of intensive studies (Rowan-Robinson 
1988, van den Bergh 1992, and de Vaucouleurs 1993). 
Broadly speaking there are two distinct methods of cali- 
brating distances to galaxies. The first group of methods 
relies on accurate distances (parallaxes) to nearby objects 
to calibrate a 'distance ladder' extending to objects further 



away. A recent successful example of this is the HST Key 
project (Freedman et al. 1994 and Freedman et al 1997) 
that aims to measure the Hubble constant with an accuracy 
of 10% by using Cepheid variables as standard candles to 
measure distances to the Virgo Cluster. Current estimates 
of Ho using this and similar methods (Tanvir et al. 1995) 
are in the range w 60-90 kms -1 Mpc" 1 . By accurately cali- 
brating distances to the Virgo cluster galaxies the 'distance 
ladder' will be extended via numerous secondary methods, 
including the D-cr relation (Faber et al. 1989), Tully- Fisher 
relationship (Tully & Fisher 1977), surface-brightness fluc- 
tuation method (Tonry & Schneider 1988), and supernova 
(Sandage et al. 1996 and Perlmutter et al. 1997) method. 

The second group of methods does not require any cal- 
ibration or progression along a 'distance ladder' but applies 
directly to the object concerned. Examples include the use 
of gravitational lens systems (Refsdael 1964 and Kundic et 
al. 1997) and the Sunyaev-Zel'dovich effect ( Sunyaev & 
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Zel'dovich 1980 and McHardy et al. 1990). Current early 
estimates of Ho using these two methods are in the range 
« 30-80 kms -1 Mpc" 1 . These direct methods, whilst gen- 
erally more model dependent, give an important check on 
'distance ladder' methods. 

AGNs are potentially important cosmological probes 
because their high luminosity allows them to be observed 
at large redshifts. Correlations between the luminosity and 
various emission line strengths and ratios (Baldwin 1977 and 
Kinney et al. 1990) have been investigated for many years, 
but have not yet allowed a consistently accurate inference 
of the distances of AGNs independently of their redshifts. 
A promising geometrical method that uses proper motions 
and line-of-sight accelerations of water maser emission in 
NGC 4258 yields an accurate distance of 6.4 ± 0.9 Mpc 
(Miyoshi et al. 1995). In principle, a similar method may 
be employed to infer distances to maser sources but not 
to much higher redshifts where proper motions become too 
small. 

We propose a new method that utilizes the relatively 
simple physics of light-travel time and blackbody radia- 
tion to measure directly redshift-independent luminosity dis- 
tances to AGNs, and hence determine Ho- In §2 we dis- 
cuss the theory of the method, and show how wavelength- 
dependent time delays determine T(R), and measure Ho- 
Therefore, wavelength-dependent time delays can be used 
to a) provide 'smoking gun' evidence of accretion disc struc- 
tures in AGN and b) determine cosmological parameters, 
e.g Ho- The first application of our method uses data 
from the International AGN Watch monitoring campaign 
on NGC 7469, and is discussed in §3. The systematic errors 
in our method are presented in §4, and §5 summarizes the 
main results of the paper. 



2 THEORY 

A blackbody accretion disc illuminated by a central source 
has a radial temperature profile T(R) that is a non-linear 
combination of the surface temperature due to viscous heat 
dissipation T v i s and that due to irradiation Ti rr of the disc, 
T 4 = T 4 S + 7] 4 r . Therefore T(R) depends on both the ge- 
ometry of the accretion disc and the relative prominence of 
viscous heat dissipation and irradiation effects. For exam- 
ple, when T(R) is determined by viscous dissipation alone, 
a T oc (MM) 1/4 R~ 3/i structure exists where M is the mass 
of the black hole and M is the mass accretion rate. An irra- 
diating source, luminosity L, situated a height // x above the 
disc plane incurs a similar T oc {LHx) l '*R~ 3 ' 4 structure for 
R 2> -ff x , provided the disc thickness H <C H x . 

The reprocessing hypothesis assumes that the 
UV /optical continuum variations represent the response of 
gaseous material to variations in the higher-energy contin- 
uum. The stringent upper limits, < 0.3 day, on time delays 
between the X-ray and UV variations in NGC 4151 (Edelson 
et al. 1996) suggest that the variations in different wave- 
bands must be radiatively coupled (i.e., any possible time 
delays are due to light-travel time effects), since, for exam- 
ple, viscous time scales are much too long. Furthermore, the 
equivalent width of Fe Ka at 6.4 KeV and the strength of 
Compton reflection observed at > 10 KeV suggest that the 
majority of X-rays generated by an isotropic source must 



be reprocessed by relatively cold (< 10 6 K) optically thick 
gas, possibly that of an accretion disc (Pounds et al. 1990 
and George & Fabian 1991). This has led to 'ad hoc' models 
where the higher-energy, e.g X-ray, source illuminates the 
accretion disc from above. 

We assume here that a similar mechanism must be op- 
crating, i.e. some variable source of high-energy radiation 
in the vicinity of the disc axis illuminates the disc and ra- 
diatively drives the UV/optical continuum variations. When 
high-energy radiation is emitted from the central regions of 
the disc, a wave of heating propagates out at a speed c arriv- 
ing at radius R after a mean time t — R/c. At this radius R 
the temperature T(R) rises slightly, thereby emitting more 
photons near wavelength A = hc/kTX (where X ~3-4 for 
blackbody radiation). When we observe a time delay r at 
wavelength A that in effect measures the radius R = re at 
which the disc has temperature T = hc/kXX. 

Assuming a temperature profile of the disc, T — 
To(-R/i?o) _3//4 , the wavelength-dependent time delay is 

r(A)= 3.9d(^) 4/3 (^) V3 (f) 4/3 , (1) 

where To is the temperature of the disc at radius Ro = 1 
light day. Note that a disc with T oc i? _3//4 predicts a r oc 
A 4 / 3 wavelength-dependent time delay, and an annulus of 
fixed R in the disc responds with a range of time delays 
r = (R/c)(l ± sini), with i the disc inclination. Hence, in 
principle, the width of the time delay distribution at each 
wavelength determines i. The explicit inclusion of X in the 
above and following equations is for heuristic purposes, since 
X is not a free parameter or variable but is determined by 
the blackbody model. 

With T(R) determined from the observed r(A), the pre- 
dicted spectrum of the disc can be calculated straightfor- 
wardly by summing up the blackbody contributions from 
various disc annuli (Shakura and Sunyaev 1973). The disc 
spectrum is given by 

U = 11.2 Jy (^) 2 (^)- 2 (^)" 3 (f )-"««, 

where D is the distance to the AGN and i is the incli- 
nation of the disc. Note that the classical thin disc spec- 
trum, /„ oc A -1 / 3 , is recovered since r oc A 4/3 , hence 
/„ oc r 2 A -3 oc A -1 / 3 . We note that /„ is the distribution 
of flux with frequency v. 

The redshift-independent distance to the AGN is then 
derived to be 

D = 3.3 Mpc (^) ^ (^r /2 (f )" 4/3 

By inserting observed values of r(A) and /„ into the above 
equation, we determine a redshift-independent distance to 
the object. Hubble's constant, Ho = cz/D, is then, 

*o= 89.6 ^ (^) 3/ \o^) i^) 1 "- 
(f) 4/3 - 

where z « 1 is the redshift of the AGN. 
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Figure 1. Top Panel: The predicted time delays for the irradiated accretion disc model are compared with the observed time delays, 
measured relative to 1315A. The time delay increases above the overall trend near the emission lines because here there is a mix of 
continuum and line flux. The lines respond with longer time delays and this results in a larger net delay. Bottom panel: Model spectra 
for an irradiated accretion disc (assuming Ho = 42kms~ 1 Mpc -1 and i = 45°) are compared with observed spectra of NGC 7469 from 
the 1996 AGN Watch multiwavelength monitoring campaign. See text for full details. 



3 APPLICATION OF METHOD TO NGC 7469 
MONITORING DATA 

To apply our new method we use data from a seven-week 
International AGN Watch (Alloin et al. 1994) multiwave- 
length monitoring campaign on NGC 7469, z = 0.0164, 
which showed optical continuum variations lagging behind 
those in the UV by about 1 day at 4800A and about 2 days at 
7500A (Wanders et al. 1997 and Collier et al. 1998). These 
UV /optical continuum lags have been shown to be statisti- 
cally significant at no less than 97% confidence (Peterson et 
al. 1998). The measured time delay between the flux vari- 
ations at wavelength A and those at wavelength 1315A is 
shown in the top panel of Figure 1. The observed time de- 
lays are encompassed by the r oc A 4 / 3 predictions, shown 
as dashed and dotted lines, for irradiated blackbody accre- 
tion disc models with T = To(R/ Ro)~ 3 ^ 4 structure, with 
T min « 6500 K and T max « 7700 K for faint and bright 
states respectively. The observed delays increase above the 
overall trend at wavelengths near emission lines. Here there 
is a mix of continuum and line flux, the lines responding 
with larger delays than the continuum light. The horizontal 
axes in Fig. 1 give for each wavelength A the corresponding 
temperature T in the disc. The temperature is calculated as- 
suming T = hc/k\X, with X = 3.89 being appropriate for 
blackbody discs with T oc J? -3 / 4 as indicated by numerical 
simulations. 

The bottom panel of Fig. 1 compares the de-reddened 



UV/optical maximum, minimum and difference spectra of 
the NGC 7469 monitoring campaign with predicted spectra 
for Ho = 42 kms" 1 Mpc" 1 (D « 117 Mpc) based on the 
irradiated accretion disc model inclined at i = 45° to the 
observer's line-of-sight. The dotted and dashed lines repre- 
sent model spectra for bright and faint states of the irra- 
diated accretion disc, and define maximum and minimum 
temperatures, T max « 7700 K and r min « 6500 K, at ra- 
dius Ro = 1 light day respectively. The brightest and faintest 
spectra seen during the NGC 7469 campaign are both much 
redder than the predicted disc spectra. This is caused by 
contamination of the observed spectra by a red starlight 
component from the host galaxy. We therefore consider these 
spectra to be upper limits to the spectrum of the active nu- 
cleus. The solid line represents the predicted difference spec- 
trum between the bright and faint states of the disc. This 
agrees approximately with the difference spectrum between 
the brightest and faintest spectra recorded in the NGC 7469 
AGN Watch campaign. The difference spectrum, which can- 
cels any starlight contamination, gives a lower limit to the 
nuclear spectrum. 

The results shown in Fig. 1 demonstrate that the ob- 
served variability in the continuum spectrum of NGC 7469 
is in approximate agreement with a blackbody disc, T oc 
R-°- 75 , and H (cosi) 1/2 = SSiekms -1 Mpc -1 . The error- 
bars reflect statistical uncertainties in the measured fluxes 
/„, time delays r and redshift z, e.g. equation (4). The point- 
to-point scatter in the continuum time delays allows us to as- 



© 1994 RAS, MNRAS 000, Q-?? 



4 S.J. Collier et al. 



sign a 10% uncertainty to the time delay measurements. The 
lower limit to the flux of the nucleus, is A/„ = 3.5± 1.0 mjy 
at A = 1315A. The uncertainty in the redshift of NGC 7469 
(z — 0.0164), induced by assuming ±300 km s" 1 peculiar 
velocities (Lynden-Bell et al. 1988), is 6%. Our 10% uncer- 
tainty in r, 29% and 6% uncertainties in A/„ and z respec- 
tively, result in a 18% uncertainty in Ho (cos i) 1 ^ 2 . 



4 SYSTEMATIC ERRORS 

Our estimate of Ho is subject to several systematic errors. 
The AGN spectrum is diminished and reddened by interven- 
ing dust. Reddening estimates derived from pointed 21cm 
observations give E(B-V) w 0.074-0.096 (Elvis et al. 1989 
and Lockman & Savage 1995). Other estimates based on 
using the 2200A dust absorption feature give E(B — V) ~ 
0.14 (Westin 1985). We have corrected our spectra using 
E(B — V) = 0.14. The host-galaxy contamination can be 
estimated from off-nuclear observations of the host galaxy, 
e.g using HST or ground based adaptive optics. In the dif- 
ference spectra the host-galaxy contamination is negligible. 
The red slope of the mean spectrum is due to contamination 
by starlight from the host galaxy, which contributes at least 
« 40% at 5400A in a circular 10 arcsecond radius aperture 
(Malkan & Filippenko 1983). Welsh et al. 1998 present con- 
temporaneous HST observations of NGC 7469 and estimate 
a percentage host galaxy contamination at 7400A of ~ 80% 
in a 10" x 16".8 aperture. 

A systematic error arises from uncertainty in (cosi) 1//2 . 
However, the inclination uncertainty is not a major obsta- 
cle. According to unified schemes (Antonucci 1993 and Hes 
et al. 1993), Seyfert 1 galaxies, in which we see the broad 
emission line region (BLR), have i < 60°, while Seyfert 2 
galaxies, in which the BLR is obscured by a dusty torus, 
have i > 60°. For i < 60°, (cos i) 1/2 > 0.7. By adopting 
i — 45°, (cosi) 1 / 2 = 0.84, we commit a maximum error 
of ± 17%, and an RMS error of ± 11%. Averaging over 10 
objects could reduce this by a factor \/T0 to ± 3.5%. With 
(cos i) 1/2 = 0.84 ± 0.1 we find H = 42 ± 9kms _1 Mpc -1 . 
We may be able to reduce the uncertainty in (cosi) 1 ^ 2 if 
we can measure i independently rather than averaging over 
the full range of i. Our echo mapping method in principle 
allows us to derive the inclination i, because the width of 
the time delay distribution at each wavelength is a function 
of i. This method may be applied, in future, to NGC 7469 
and other Seyfert 1 galaxies. However, it is likely to be a 
difficult task. Fits to the profile of the X-ray FeKa line in 
MCG-6-30-15 (Tanaka et al. 1995) yielded a disc inclina- 
tion i of 30 ± 3°. Similar observations of an X-ray FeKa 
line in NGC 7469 could therefore measure i with about 10% 
accuracy. It may also be possible to determine i from polar- 
ization measurements. Therefore there are good prospects 
for measuring inclinations of individual AGNs. 

The blackbody model is a source of systematic uncer- 
tainty. Our value X — 3.89 relies on the assumption that 
the changes in the UV/optical continuum can be modelled 
as irradiation of a blackbody disc. This is justified by the ap- 
proximate agreement of the predicted /„ oc A" 1//3 spectrum 
with the observed difference spectrum. This can be further 
investigated by considering models of the vertical structure 
and the emitted spectra of irradiated accretion disc atmo- 



spheres (Sincell & Krolik 1997) . We expect that to first order 
our method of measuring distances should be insensitive to 
limb darkening (Hubeny et al. 1997). First, the irradiation 
of the disc will flatten the temperature versus optical depth 
relationship in the atmosphere and conspire to reduce the 
limb darkening effect. Second, the lower temperatures ob- 
served at high inclinations will change the apparent T(R) 
profile, but the same blackbody relationship will describe 
the surface brightness distribution. Hence, while our T(R) 
profile is sensitive to limb darkening the inferred surface 
brightness and subsequent distance is not. 

Another source of systematic uncertainty is the source 
geometry and nature of the continuum variations. The repro- 
cessing geometry we have considered might be completely 
ruled out by the recent X-ray observations of NGC 7469 
(Nandra et al. 1998). The 2-10 KeV X-ray variations are 
poorly correlated or uncorrelated with the UV/optical vari- 
ations described here. Either (a) the accretion disc sees dif- 
ferent X-ray variations than we do, (b) the UV/optical vari- 
ations are driven by another unobserved part of the spec- 
trum, e.g. the extreme UV, or (c) the model geometry is 
completely wrong. However, observations of time delays be- 
tween different wavebands and knowledge of the spectral 
energy distribution of sources will allow us to determine the 
importance or irrelevance of these various systematic effects. 



5 SUMMARY 

Our Hubble constant estimate of Ho = 42±9kms~ 1 Mpc" 
is not consistent with the independent 'distance ladder' es- 
timates of Ho — 80 ± 17kms _1 Mpc -1 by Freedman et al. 
1994, and H = 69 ± 8kms _1 Mpc" 1 by Tanvir et al. 1995. 
However, the 22% accuracy of our Hubble constant estimate 
compares favourably with the 21% and 12% uncertainties of 
the Hubble constant estimates reported by Freedman and 
Tanvir respectively. For comparison, the direct method uti- 
lizing gravitational lenses give, for example, Hubble con- 
stant estimates of Ho = 42 ± 6kms~ 1 Mpc" 1 (Schecter et 
al. 1997), Ho = 53^7° kms" 1 Mpc" 1 (Courbin et al. 1997), 
and H = 64 ± 13 kms" 1 Mpc" 1 (Kundic et al. 1997). The 
Sunyaev-Zel'dovich method gives Ho = 47l 2 jkms- 1 Mp C - 1 
(Hughes and Birkinshaw 1998 and references therein). Any 
apparent discrepancy in the Hubble constant estimates from 
'distance ladder' and direct methods will need investigation, 
which in turn requires statistically significant samples of the 
methods to be compared. 

Our result based on observations not specifically de- 
signed to measure Ho can be improved upon. A continu- 
ous 2-3 month multiwavelength, multi-telescope monitoring 
campaign on a sample of AGNs will make this method a 
serious competitor with the established methods of measur- 
ing Ho- Principally, the time delay measurements can be 
constrained to better than 5%. There are good prospects 
for measuring inclinations of individual AGNs with high ac- 
curacy. The starlight contamination can be estimated and 
corrected for as already described, allowing an accurate de- 
termination of the nuclear spectrum to better than 10%. 
Finally, a statistical average of individual Ho measurements 
will reduce our final uncertainty, and we expect a future op- 
timally designed experiment to determine the Hubble con- 
stant with < 10% accuracy. 
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To conclude, we identify the variable component of 
the de-reddened UV and optical continuum fluxes as a 
lower limit to the nuclear spectrum. For NGC 7469 we 
note that this agrees approximately with the f v oc A -1 / 3 
spectrum predicted for a blackbody accretion disc with a 
T oc R~ 3 / A structure. At the same time the wavelength- 
dependence of the observed time delays, r oc A 4,/3 , is con- 
sistent with T oc i? _3//4 . The concurrence of these two in- 
dependent lines of evidence strongly supports the notion 
of a standard blackbody accretion disc in NGC 7469, and 
strengthens the evidence (Shields 1978, Malkan et al. 1982, 
and Tanaka et al. 1995) for accretion discs in AGN. Us- 
ing the variable component of the continuum fluxes we find 
H (cosi/0.7) 1/2 = 42±9kms _1 Mpc" 1 . The interpretation 
of the observed time delays and spectral energy distribution 
in the context of an accretion disc structure requires fur- 
ther validation. However, analysis of the observed variable 
spectrum and wavelength-dependent time delays along the 
lines outlined above yields redshift- independent luminosity 
distances to AGNs. This opens up a new route to Ho and 
by extension to fainter objects at z ~ 1, qo- 
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